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ABSTRACT 

We study how the properties of transient sources of ultra-high-energy cosmic rays (UHECRs) can 
be accessed by exploiting UHECR experiments, taking into account the propagation of UHECRs 
in magnetic structures which the sources are embedded in, i.e., clusters of galaxies and filamentary 
structures. Adopting simplified analytical models, we demonstrate that the structured extragalactic 
magnetic fields (EGMFs) play crucial roles in unveiling the properties of the transient sources. These 
EGMFs unavoidably cause significant delay in the arrival time of UHECRs as well as the Galactic 
magnetic field, even if the strength of magnetic fields in voids is zero. Then, we show that, given 
good knowledge on the structured EGMFs, UHECR observations with high statistics above 10 20 eV 
allow us to constrain the generation rate of transient UHECR sources and their energy input per 
burst, which can be compared with the rates and energy release of known astrophysical phenom- 
ena. We also demonstrate that identifying the energy dependence of the apparent number density of 
UHECR sources at the highest energies is crucial as a clue to such transient sources. Future UHECR 
experiments with extremely large exposure are required to reveal the nature of transient UHECR 
sources. 

Subject headings: cosmic rays — magnetic fields — methods: numerical 



1. INTRODUCTION 

The origin of ultra-high-energy cosmic rays (UHE- 
CRs) has been a mystery for more than forty years. 
The highest energy cosmic rays (> 10 19 eV) are usu- 
ally thought to be extragalactic origin, and various 
kinds of astrophysical objects have been suggested 
as primary source can didates, i ncludi n g gamma-ray 
bursts (GRB s ) (e.g . iWaxmanl fl995t iVietril 
IMurase et all 120061: IMurase et all I2008af). 
born magnet ars (|Aronsl l2003t IMurase et al.l 
Koteral 120111). active g al actic nu clei (AGN ) 
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shocks (e.g.. iNorman et al.l 119951: iKang et al.l 119961: 
llnoue et all 12007ft . Theoretically. UHECR sources are 
expected to be powerful enough. For cosmic-ray accel- 
erate^ associated with an outflow, the Hillas condition 
(|Hillaslll984ft can be r ewritten in terms of the iso tropic 
luminosity L as (e.g.. [Blandford 1 120001: IWaxmanl I2004t 
IFarrar fe Gruzinovll2009l; ILemoine fe Waxmanll2009D 



L B = e B L > 2 x 10 45 



Z 2 f3 



erg s 



(1) 



where e B , Z, T, j3 and E 2 q = E/10 eV are a frac- 
tion of magnetic luminosity to the total luminosity, the 
nuclear mass number of cosmic rays, the bulk Lorentz 
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factor of the outflow, the velocity of a shock or wave in 
the production region in the unit of speed of light and 
the energy of cosmic rays, respectively. Among known 
candidates, few steady sources such as Fanaroff-Riley 
(FR) II galaxies seem to satisfy this condition in local 
Universe for Z = 1, which is inconsistent with the ob- 
served anisotropy as long as UHECRs are protons (e.g., 
iTakami fe SatoTl2009f ) . Also, Zaw et al. (2009) argued 
that the power of AGN correlating with detected UHE- 
CRs seems insufficient to produce UHECR protons. The 
above luminosity requirement can be satisfied, however, 
if UHECRs are generated by powerful transient phenom- 
ena like AGN flares, GRBs a nd newly born magnetars 
even if they are protons (e.g.. IFarrar fc Gruzinovl [20091 : 
iDermer et aT1l2009l: ILemoine fc Waxmanll2009f )~ 

The other possible astrophysical solution is to con- 
sider that heavy nuclei dominate over protons, where 
the required luminosity is reduced by Z 2 and therefore 
more objects are allowed to be UHECR sources. In- 
deed, the heavy-ion-dominated composition has been im- 
plied b y recent results of the Pierre Auger Observatory 
(PAO) (| Abraham et al.ll2010al) . If this is the case, only a 
few nearby radio galaxies or even a single AGN such as 
Cen A may contribute to the observed UHECR flux (e.g., 
iGorbunov et al.ll2008f). Other sources, including radio- 



quiet AGN dPe'er et al1 l200l and GRBs (|Murase et al.l 
120081 IWang et alJ l2008) , are also viable. The absence of 



anisotrop y at ~ 10 20 eV/Z m ay imply high abundance 
of nuclei (Lemoinc fc Waxmari 12001 lAbreu et all 12011 
even at thelower energies, the origin of which is unclear. 
On the other hand, the PAO data on the fluctuation 
of AT max seem difficult to be reconciled with the X m ^ 
distribution of the same data (|Anchordoaui et aLl feOll). 
and proton composition may be pos sible with a different 
estim ator of primary composition (Wil k fc Wlodarczvkl 
120111 ). Also, the High Resolution Fly's Eye (HiRes) 
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has clai med proton-dominat ed composition even above 
10 19 eV (jAbbasi et al.|[20Toh . There are different argu- 
ments and the UHECR composition has not been settled 
experimentally. Proton composition seems possible at 
present. 

If UHECR sources are transient, that is, the source 
activity is shorter than the dispersion of the arrival time 
produced by cosmic magnetic fields during propagation, 
the direct identification of UHECR sources by UHECR 
observations is a more difficult task than that for steady 
sources due to the delay of the arrival time between 
UHECRs and other neutral messengers (photons, neutri- 
nos and gravitational waves) emitted by the same source 
activity. Multi-messenger approaches are definitely pow- 
erful, but it is important to extract as much information 
as possible from UHECR observations as one of the 
messengers. For transient sources, relations between ob- 
served quantities and the properties of UHECR sources 
have been discussed , considering intervening cosmic 
magnetic fields (e.g.. iMiralda-Escude fe Waxmanl [19961 : 
Wax man fc Mi ralda-Es cudel 119961 : iMurase fe Takam 



2009). Any candidate of primary UHECR sources 



has to possess enough energy budget to reproduce the 
observed flux, which is the product of the energy input 
per activity and the rate of bursts or flares. The rate 
of bursts or flares is related to the apparent UHECR 
source number density, which can be determined from 
anisotropy in UHECR arrival distribution. The relation 
between the rate and the apparent source number 
density is less obvious, depending on the Galactic 
magnetic field (GMF) and poorly known extragalactic 
magnetic fields (EGMFs). The GMF unavoidably 
affects UHECRs arriving at the Earth. Considering the 
GMF allows us to evaluate the above relation and to 
give a constraint on candidates of transient UHECR 
sources through comparing between the rate inferred 
from UHEC R observations and that of known transient 
phenomena (|Murase fc Takamil[2009l) . 

Despite current uncertainty in EGMFs, it has been be- 
lieved that the inhomogeneity of the EGMFs in the Uni- 
verse is c r ucial f or the propagation of UHECRs (e.g., 
Siel et al.1 l200l l200l iTakami et al.1 I200& iDas et alJ 



20081: iKotera fc Lemoind I2008D . The Universe indeed 



has structures, which consist of clusters of galaxies, fila- 
ments, sheets and voids. It has been suggested that mag- 
netic fields in the structured regions were amplified via 
cosmological structure formation, and various numerical 
simulations of the structure formation have shown that 
the EGMF distribut i on fo l lows the matter distribution 
(e-g JSigl et"ani200H l200i IDolag et al.ll2005t IRvu et al.1 
2008). So, astrophysical objects including astrophysical 
UHECR sources are generally embedded in the struc- 
tured regions. Thus, these magnetic structures also un- 
avoidably affect the propagation of UHECRs. The struc- 
tured EGMFs play essential roles in the time-delay and 
time-profile spread of UHECRs as well as their deflec- 
tions. 

In this paper, we study the roles of structured EGMFs 
on the propagation of UHECRs, and see how the prop- 
erties of transient UHECR sources, e.g., the rate of 
UHECR bursts or flares p s , can be constrained. In sec- 
tion^ we describe basic relations between the UHECR 
burst rate and observational quantities that can be ob- 
tained by UHECR experiments. Possible signatures of 



transient UHECR sources that may be seen in the ar- 
rival distribution of UHECRs are also discussed. In sec- 
tion[3l we numerically calculate the propagation of UHE- 
CRs in the magnetized structures for various source posi- 
tions, and evaluate the time spread due to the structured 
fields. Then, in section|4j we discuss possible constraints 
on the UHECR burst rate and cosmic-ray energy input 
burst, taking into account large uncertainty in the void 
EGMF. We present discussions in section [5] and summa- 
rize this study in section [6] Throughout the paper, the 
proton-composition is assumed, and ACDM cosmology 
with H a = 71 km s _1 Mpc -1 , Q m = 0.3 and fi A = 0.7 is 
adopted. 

2. BASIC RELATIONS 
2.1. Maximum distance of contributing sources 

UHECRs above 5 x 10 19 eV cannot avoid energy- 
loss due to photomeson production wi th cosm i c mi- 
crowave background (CM B) photons (jGreisenl 19GG; 
iZatsepin fc Kuz'minl I1966D . Hence, nearby UHECR 
sources within the so-called Greisen-Zatsepin-Kuz'min 
(GZK) radius mainly contribute to the observed flux, 
which is typically the energy-loss length of ~ 100 Mpc 
at 6 x 10 19 eV. 

For practical purpose, we here introduce D m:ix (E) that 
is the maximum distance of sources of UHECRs observed 
with the energy E at the Earth. This is defined as the dis- 
tance within which sources are responsible for 99 % of the 
total observed flux. T his definition is mo r e app ropriate 
than what we used in IMurase fc Takami (200y), where 
the energy-loss length of the photomeson production is 
simply adopted, because sources outside the energy-loss 
length can still significantly contribute to the total ob- 
served flux due to a volume effect. We estimate D max (E) 
through a backtracking method of the propagation of 
UHECRs with a continuous energy-loss approximation 
on the photomeson producti on and Bethe-Heitl er pair 
creation with CMB photons (jTakami et al.1 12006). Fig- 
ured] shows D max (E) calculated on the assmption that 

UHECR emissivity is proportional to (1 -I- z) 3 and the 
spectral index of -2.6, where z is cosmological redshift. 
Note that D max (E) is not sensitive on these assumptions 
since it is essentially determined by the energy-loss of 
the photomeson production. One sees that D max (E) 
rapidly decreases owing to the photomeson production 
with CMB photons with increasing the energy of protons. 
Anax(£) is typically 200 Mpc and 75 Mpc at E = 6x 10 19 



eV and 10 20 eV for E, 



10 eV, respectively. 



2.2. Deflection and time profile spread 

Trajectories of UHECRs are deflected in magnetized 
extragalactic space, whose m agnetic fields are still un- 
certain (e.g., iKronbera 119941 for a review). Assuming 
that EGMFs are incoherent (A <C D), a typical deflec- 
tion angle of UHECRs durin g propagation in extragalac- 
tic sp ace is estimated to be (jWaxman fc Miralda-Escudel 
fl9M ) 

-2° B_ 9 \ 1/2 E 20 - 1 D 75 1 / 2 
<20°E 2Q ^ 1 D 75 1 / 2 , (2) 
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Figure 1. Maximum distance of UHECR sources which can con- 
tribute to the flux observed at the Earth at the observed proton 
energy E on the assumption of -E max = 10 21 eV (solid) and 10 20 ' 5 
eV (dotted). The maximum distance rapidly decreases due to pho- 
tomeson production with CMB photons with increase of the proton 
energy. Note that £> ma x(E) is 200 Mpc at E = 6 X 10 19 eV and 
75 Mpc at E = 10 20 eV for £ max = 10 21 eV. 



where -D75 = D/75 Mpc and n, = E/eB are the propa- 
gation distance of UHECRs and the Larmor radius of the 
UHECRs for the characteristic EGMF strength of B_ 9 = 
B/10- 9 G and the correlation length of A = A/10° Mpc, 
respectively, e is the electron charge magnitude. The 
last inequality is obtained from an upper limit on the 
averaged EGMF from Faraday rotation measurements of 
distant quasars, BX 1 ' 2 < (10nG)(lMpc) 1 / 2 (|Rvu et al.l 
1998; Bla si et all [1999). Similar values have been de- 
rived from CMB and matter power spe ctra for the pri- 
mordial magnetic field at p resent (e.g., Uedamzik et al.1 
120001 lYamazaki et al.l l2~010). On the other hand, several 
plausible GMF models predict deflection angles of < 5° 
at ~ 6 x 10 19 eV and smaller at higher energies in al- 
most all the direction s in the sky ()Takami fc Satol [20081 : 
Takami fc Satoll2010T l. Those deflections seem compara- 
ble to the the correlation angle between UHECRs and 
their source candidates which were ori ginally reported 
by th e PAO, - 3° above - 6 x 10 19 eV ((Abraham et al.l 
2007). However, one should keep in mind that such angu- 
lar scale should not always be interpreted as the typical 
deflection angle of UHECRs, i.e., the source candidates 
could be tracers of true sources, both of which are em- 
bedded in the large-scale structure which are spread with 
typical angular scale of ~ 10° in local Universe in the 
projected sky. A careful correlation analysis on the PAO 
data with galaxies in local Universe suggested the corre- 
lation with the large scale structure, imply ing the deflec- 
tion angles of < 15° more conservatively (1 Takami et al.l 
2009). 

UHECRs propagate along different paths in magne- 
tized space, which produce the time profile spread of a 
UHECR burst. The spread is comparable w ith the time- 
delay t d (jMiralda-Escude fc Waxmanll 19961 ). The appar- 
ent duration of a UHECR burst is 



a{E,D)~t d 



D9 2 {E,D) 
4c ' 



(3) 



This expression depends on D. Using D max (E), let 
us introduce the characteristic time spread by (e.g., 



iMurase fc Takamil [20091 ) 

r(E) = 



A»ax(i3) 



^D ma J{E) Jo 

\a{E,D m&x {E)). 
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dD 4ttD 2 <j(E, D) 



(4) 



In the above discussion, we have not specified EGMFs. 
The EGMFs consist of magnetic fields in structured re- 
gions, i.e., clusters of galaxies and filaments, and mag- 
netic fields in voids. Clusters of galaxies typically pos- 
sess ~ 0.1 — 1 //G, whereas magnetic fields in fila- 
ments are expected to be weaker. Recent sophisti- 
cated numerical simulations suggested ~ 10 nG and 
~ 30 nG as volume-av eraged and root-m ean-square com- 
ponents, respectively (jRvu et alj|2008t) . The EGMF in 
a void is much less understood. Upper limits on the 
effective EGMF have been obtained from the Faraday 
rotation measure as mentioned above, i.e., BX 1 ^ 2 < 
10~ 8 G Mpc 1 / 2 . A recent evaluation of the upper 
limit of the primordial magnetic field, which is believed 
to be the EGMF in voids, derived B v < 2.5 nG for 
A v = 1 Mpc from the pow er spectra of CMB and mat- 
ter (jYamazaki et al.l I2010T ) . On the other hand, the us- 
age of the pair- halo /echo emission from TeV blazars ma y 
give lower bounds (jPlaga 1 119951 : IMurase et al. I l2008bf ) , 
and B v > 10~ 18 — 10~ 1 7 G have recently been suggested 
from Fermi data (e.g., IDolag et al.l 120111 : iDermer et al.l 
120111: ITakahashi et alJl201lD . If the strength of the void 
EGMF is close to these lower limits, only EGMFs in 
the structure regions play the role in the propagation of 
UHECRs. The effective magnetic field strength B e g and 
its correlation length A e ff can be estimated from, 



'(E,D)~f Q 



2DX C 

9r L , c 2 



+ /f 
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9r L , f 2 



+ /v 



2L>A V 
9rW 



2DX, 



cir 



9r L , eff 2: 



( 5 ) 

where / x , A x and rL, x for x = c, f, v are the volume 
filling fraction of a structure x, the correlation length of 
a magnetic field and the Larmor radius of UHECRs in 
a structure x, respectively, and x = c, f, v correspond 
to clusters of galaxies, filamentary structures and voids, 
respectively. 

Here, let us assume that the EGMF in voids is neg- 
ligibly weak for the UHECR propagation, focusing on 
effects of the structured EGMFs. Because of the small 
volume filling fraction of galaxy clusters (/ c ~ 10~ 4 ) in 
our mo del (see the next section) and seve ral simulation 
results (IDolag et al J 120051: iDas et al.l [20081 see also Fig- 
ure 9 of iKotera fc Olintol (|201lD ). one may approximate 
that only UHECRs generated in clusters of galaxies are 
affected by a magnetic field of clus ters. Assuming h ~ 
0.01, B t = 10 nG and A f = 100 kpc (|Rvu et al.ll2008h . an 
estimated value is B c s Xcg 1 ^ 2 ~ 0.3 nG Mpc 1 / 2 because 
only the second term contributes if all UHECR sources 

1 /2 

are embedded in filamentary structures, or B c g X c g 1 ~ 
1.0 nG Mpc 1 / 2 if all UHECR sources are embedded in 
clusters of galaxies on the assumptions of B c — 0.3/iG 
and A c = 100 kpc. 

In this section, for simplicity, the energy- loss of protons 
was neglected except for D max (E) in discussing the ana- 
lytical expressions. The energy-loss during propagation 
will be taken into account when numerical calculations 
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are performed in section [3] and 21 

2.3. UHECR burst rate and apparent source density 

We call sources "transients" when the intrinsic dura- 
tion of UHECR production at a source ST is shorter than 
the characteristic time profile spread t(E) at the ob- 
served energy of E. If the time dispersion is longer than 
the time scale of UHECR observations, we misperceive 
that UHECR bursts are steady sources, and therefore can 
define the " apparent" number density of UHECR sources 
n s (E). The source number density is related to the rate 
of U HECR bursts p s as (iMiralda-Escude fc Waxmanl 
fl99fih 

The source number density generally depends on 
UHECR energies, since the apparent duration is depen- 
dent on energies explicitly. 

UHECRs observed at the Earth suffer from the GMF 
and EGMFs embedding their sources. The GMF typically 
has order of pG for a disk component, and it also has 
random and halo components. In principle, p a can be es- 
timated by equation ([5]) if the time spread by these fields 
and the EGMF in voids can be well estimated. However, 
the EGMF in voids is highly uncertain as discussed in 
the last subsection, but could contribute to the total time 
spread significantly because of large propagation distance 
compared to the size of Galactic space and the magnetic 
structures around sources. This uncertainty leads to a 
finite range of allowed values of p s . Given inevitable con- 
tributions of the GMF and EGMFs embedding UHECR 
sources to the apparent duration, r m i n (E), and the al- 
lowed maximal time spread including the contribution 
from the poorly known EGMF in v oids, T ma ^(E), the 
rate of UHECR bursts is limited as (|Murase fc Takarml 
2009) 



n s (E) n s {E) 



(EY 



(7) 



Here, n s (E) can be, in principle, estimated from 
anisotropy in the arriva l dis t ribution of UHE- 
CRs (e.g.. [Yoshiguchi et all [200l iTakami et all [20061: 



ITakami & Sato 2009; Cuoco et all 12009). assuming that 
the time spread is longer than the UHECR observation 
timescale. 

However, one should keep in mind that equation (|6]) is 
valid when each UHECR burst can be ind ividually iden- 
tified as a burst (|Murase fc Takami 2009). If more than 
one bursts or flares occurring in an angular patch con- 
tribute to UHECRs observed in the same time-window, 
i.e., the time profiles of two independent UHECR bursts 
from the same direction (within the size of the angular 
patch) are overlapped at the Earth, equation ^ cannot 
be used as it is. Therefore, one has to focus on UHE- 
CRs with higher energies to examine cases where t(E) is 
shorter than the apparent time interval between bursts 
or flares occurring in the same angular patch, AT. In 
reality, UHECRs have finite deviation angles due to cos- 
mic magnetic fields, so UHECRs from a source arrives 
within a finite solid angle Ail = iti/j 2 around the source, 
which can be regarded as the appropriate size of the fi- 
nite angular patch. For a given p s , the apparent time 



interval between bursts in the region of the sky with Ail 
is estimated to be 



AT' 



AfiAnax(£)> a 



i / D max (E) 



-3 



~ 3 x 10 V s -V»-^=j y, (8) 

where -05 = V'/S and p s , — p s /!0° Gpc~ 3 yr _1 . We 
take the typical positional correlation scale as ip, and 
use ip ~ 5° as a reference choice, which corresponds to 
B eS X cS 1/2 < 2 nG Mpc 1 / 2 . This is reasonable, since 
this is consistent with the effective EGMFs estimated 
in the last subsection and a current upper limit of the 
void EGMF from a pl ausible cosmological m odel is 2.5 
nG for A v = 1 Mpc (lYamazak i et al.ll2010f ). but more 
conservative discussions with larger values of ip are also 
possible. Equation (JSJ implies that a smaller Ail gives 
larger AT, but AT should be limited to the burst/flare 
intermittence in a host galaxy, ps n^j p s , where rih is the 
number density of host galaxies of UHECR sources. In 
other words, Ail smaller than the corresponding lower 
limit is meaningless, at which one host galaxy should 
exist in a volume with a solid angle Ail within D ma , K (E). 

We call the case "bursting case" that only one burst 
or flare contributes to arriving UHECRs at a time in 
a direction, i.e., t(E) < AT. Then, the requirement 
t(E) < AT gives a sufficient condition to apply equa- 
tion ([6]), which leads to 



,(£) < 3 x HTtys" 



75Mpc 



Mpc" 



(9) 



with the usage of equation ^ . As demonstrated in Fig- 
ure [5J the range of n s (E), in which equation ^ can 
be applied, is extended to larger n s (E) at higher ener- 
gies, because the smaller number of sources contributing 
to the observed flux decreases the probability that two 
bursts are temporally overlapped in a region of the sky 
(see also equation ©). Thus, we especially focus on 
cases of E ~ 10 20 eV to demonstrate constraints on p s 
in Sections [3] and [4j although discussions are general for 
other E. Note that, although we here fix ip even at higher 
energies, smaller values of ip are expected there, so that 
the extension of the curve to higher energies would be 
more easily justified. 

On the other hand, if r{E) > AT, another UHECR 
burst may start to contribute before the end of the former 
UHECR burst is observed, and equation ([8]) implies, 



/) s >3x 10 



DmaxjE) 

t(E) V 75Mpc 



5^5" 



Gpc 



3 yr- 1 . 



(10) 



Since n s (E) can be determined by the auto-correlation 
analysis, equation © enables us to estimate p s from 
observational quantities, if t(E) can be evaluated by 
EGMF simulations and observations. Importantly, 
n s (E) has the characteristic energy dependence, which is 
demonstrated in Figure [5] Here, the case that only the 
EGMF in filamentary structures affects the time spread 
of UHECR bursts is considered for demonstration, i.e., 

-BeffAeff 1 ^ 2 ~ 0.3 nG Mpc -3 . Two representative cases 
for p s are shown, i.e., p s = 1 Gpc -3 yr -1 and 10 2 Gpc -3 
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Figure 2. The diagram of how transient UHECR sources are ob- 
served by UHECR experiments. In the region below the upper solid 
line each UHECR burst is observed with spatial and temporal sepa- 
ration, i.e., equation J5} is valid. This case is called "bursting case" 
in this paper. Note that ip = 5° is assumed. In the region below 
the lower solid line, there is no UHECR source within -D max (.E). 
The relation between the apparent source number density n B (E) 
and E is demonstrated (dashed lines) for p s = 1 Gpc -3 yr — 1 
and p a = 10 2 Gpc -3 yr _1 in the case where a EGMF in fila- 
mentary structures dominantly affects the propagation of UHE- 
CRs, i.e., B eH \ ff 1 / 2 ~ 0.3 nG Mpc" 3 . The curve in the case of 
p a = 10 2 Gpc -3 yr — 1 is saturated at which is assumed to be 
f0" 4 Mpc A 

yr _1 . One sees that n s (E) changes by more than one 
order of magnitude when E increases by the cubic root 
of ten. This means that anisotropy features are different 
among energies. Thus, observations of UHECRs above 
10 20 eV are crucial to identify this tendency clearly. Fu- 
ture UHECR experiments with large exposures may de- 
tect a large number of the highest-energy events and al- 
low us to determine the dependence of n s (E). 

Now, remember that, even if ip was small, AT would 
be limited by rih/ p s - In other words, n s (E) must not 
be larger than n^, and thus the curve of n s (E) may be 
saturated. This situation is also demonstrated, assuming 
rih = 10" 4 Mpc -3 , which is c omparable to the local n um- 
ber density of FR I galaxies (jPadovani fc Urrvlll990D . In 
the case of p s = 10 2 Gpc -3 yr _1 , the curve of n s (E) is 
saturated and then becomes flat at low energies. Thus, 
the number density of host galaxies could also be esti- 
mated if such a saturated curve is seen. Note that rih 
does not depend on a given EGMF model, and gives a 
robust constraint on UHECR source population through 
relations between UHECR sources and th eir host galax- 
ies. T he relations were briefly discussed in lTakami et al.l 

(Emi). 

There is a region on the diagram, where less than one 
transient source contributes to the total flux observed at 
the Earth in the bursting case, especially at high energies. 
This is obviously expected from the fact that D max (E) 
becomes smaller at higher energies. This region cannot 
be used to constrain UHECR sources because there is no 
UHECR source in the sky. The condition that less than 
one transient source contributes to the total flux is 

4-7T q 

—D max 3 (E)p s a(E, D max (E)) < 1. (11) 
Remembering n s (E) ~ 3p s a(E, D max (E))/5, the region 



where less than one source exists within D max (E) is de- 
termined by 

n , (E) <3xl0-'(%^grf Mpc-. ,12, 

This border line is also shown in Figure [3J below which 
the curves of n s (E) are truncated. This region does not 
depend on the choice of tp. Just above the border line, 
there are only a few sources within D max (E), where we 
can observe strong anisotropy from these nearby rare 
sources, which is strong evidence for the location of 
UHECR sources. In this viewpoint large UHECR ex- 
periments for UHECRs above E > 10 20 eV are required. 

The energy dependence of n s (E) may originate from 
transient source scenarios, but it is not always direct ev- 
idence on transient UHECR sources, although the de- 
pendence is generally different between transient scenar- 
ios and steady scenarios. If all the UHECR sources 
are identical and steady, n s (E) is a constant up to 
E = E max . However, if steady sources have different 
E max a mong sources, n^E) could de pend on E. For in- 
stance, iKachelriess fe Semikozl (|2006l ) demonstrated that 
steady UHECR sources with the power-law spectral slope 
of ~ 2.0 can reproduce the observed steep spectrum 
if a power-law distribution of E max with the index of 
~ 1.7 among sources is assumed. In this case n s (E) 
is expected to be proportional to E . On the other 
hand, in the cases of transient sources demonstrated 
in Figure [2j the spectral index of the dependence of 
n s (E) on E is much less than —2.0 above 10 20 eV ow- 
ing to the dependence of D max (E) on E. Different E max 
among sources is expected to lead to additional steep- 
ening of the dependence, although an identical E max is 
assumed throughout this paper for simplicity. There- 
fore, transient cases has the steeper dependence on E, 
but careful discussions are required to distinguish tran- 
sient cases from steady cases. Spiky spectral features 
in the spectrum of UHECRs which are specific in tran- 
sient sources may also help suppo rt transient scenarios 
(|Miralda-Escude fc Waxmanlll996[ ). 

We have shown that the dependence of n s (E) on 
E provides us with a hint of transient sources. Al- 
though above discussions neglect the energy loss of 
UHECRs, the tendency that n s (E) decreases at high 
energies is a general feature. In order to esti- 
mate n s (E) at sufficiently high energies, large statis- 
tics beyond the GZK energy (~ 10 20 eV) is re- 
quired. Large statistics can also red uce the uncer- 
tainty of estimated n s (E) significantly (|Takami fc Satol 
|2008f ). Future large UHECR experiments to target 
UHECRs above 10 20 eV such as Extreme Universe 
Space Observatory (JEM-EUSO) (jEbisuzaki et al.ll2008l) 
and th e Northern site of the Pierre Auger Obser- 
vatory (jBlumer fc the Pierre Aug er Collabor ationll20Tof) 
will give important information to constrain the proper- 
ties of transient UHECR sources. 

3. PROPAGATION OF UHECRS IN MAGNETIZED 
REGIONS 

The positions of UHECR sources in structured regions 
depend on source candidates. Recent studies have shown 
that the star formation rate of galaxies tends to be larger 
in the outskirts of clusters and filaments than that inside 
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Figure 3. Number densities of cxtragalactic backgr ound light in - 
cluding the CMB (the low infrared model of Knciskc et al.l H2004 P 
(red) and of background light in a cluster of galaxies at 100 kpc 
(green), 300 kpc (blue), 500 kpc (magenta), 1 Mpc (light blue), and 
2 Mpc (black) away from the center of the cluster. 



clusters (e.g.. iPorter et al.l feOOS). Since the explosion of 
massive stars is expected to happen in galaxies with high 
star formation rates, UHECR source candidates such as 
GRBs and newly born magnetars are more likely to be 
located in filaments or outskirts of clusters. On the other 
hand, powerful radio galaxies may be more centrally dis- 
tributed (e.g.. ILin fc Mohrl [200% . Here, two EGMFs, 
i.e., the EGMF in a filamentary structure and the EGMF 
in a cluster of galaxies, and three locations of UHECR 
sources are investigated: in a filamentary structure, at 
the center of a cluster of galaxies and 1 Mpc away from 
the center of a cluster of galaxies. In addition, the GMF 
is also considered as an unavoidable magnetic field for 
UHECRs arriving at the Earth. 

3.1. Models and calculation method 

Here, a cluster of galaxies is modeled as a spherical 
structure with the radius of 3 Mpc. Since the local num- 
ber d ensity of clusters is ~ 10~ 6 Mpc -3 (|Mazure et al.l 
1996), the volume filling fraction of galaxy clusters in this 
model is ~ 10~ 4 , which is the referred value in section 



The strength of a magnetic field in the cluster is as- 
sumed to scale with the distance from the center of the 
cluster as given by the flux-freezing condition with an 
thermal electron compo nent modeled a s a /3- model fol- 
lowing the treatment of Ide Marco et ail (|2006l ). 



B(r) =B (l 



-0.7 



(13) 



where r c — 378 kpc and B = 1/iG. The direction of the 
magnetic field is set to be turbulent with the Kolmogorov 
power spectrum with the maximum scale of A Cjmax . The 
strength of the magnetic field averaged over this cluster 
model is ~ 0.3/LtG. 

An infrared background photon field in the cluster of 
galaxies is modeled as the superposition of the spectral 
energy distributio n of 100 giant elli ptical galaxies calcu- 
lated by GRASIL (jSilva et al.lll998l) in addition to a low 
inf rared model of extragalactic background light at z — 
bv lKneiske et al.l ([2004) . The galaxies are assumed to be 



distributed follo wing an analytical fit ting formula of the 
gas distribution (|Rordorf et al.l 120041 ) . 



f(r) oc 









«1 






OL2 ' 




1 + 












oc 






1 + 







1 



(14) 

where r\ — 10 kpc, r 2 = 250 kpc, r 3 — 1 Mpc, ot\ — 
—0.51, a2 = —0.72, and a 3 = —0.58. Figure [3] shows the 
number densities of photons in the cluster of galaxies at 
several radius from the center. 

As mentioned above, the two positions of a UHECR 
source in the cluster of galaxies are considered. In both 
cases, UHECRs are injected as a jet in the direction of 
an observer, i.e., in the radial direction, with the opening 
angle of 0.1 radian (~ 6°). The jet-like UHECR injection 
is motivated by many source candidates such as GRBs 
and AGN. 

A filamentary structure is approximated to be a cylin- 
der with the radius of 2 Mpc with the magnetic strengt h 
of 10 nG based on simulation results (|Rvu et al.l [2008). 
Although the direction of a magnetic field is not clear yet, 
we assume a turbulent field with the Kolmogorov spec- 
trum with the maximum scale of Af jlnax . The photon 
field in the filament is assumed to be the same model of 
extragalactic background light as that used for the clus- 
ter. A source is located on the axis of the cylindrical 
filament and emits UHECRs toward a direction perpen- 
dicular to the axis with the jet opening angle of 0.1 ra- 
dian. Although the direction of the jet with regard to 
the axis of the cylindrical filament depends on individual 
sources, this configuration is justified because it provides 
with a reasonable lower limit of the time-delay and resul- 
tant time spread of UHECRs. In addition, this allows us 
to avoid the dependence of results on the height of the 
cylinder, although it is assumed to be 25 Mpc. 

Regarding the correlation length of the magnetic fields, 
the energy of turbulence is, in general, injected at large 
scale characterized by the driving scale of a turbulent 
magnetic field. The energy is transferred into smaller 
scale turbulence, and then is dissipated to thermal en- 
ergy at small scale. In the intermediate scale, the energy 
spectrum of the turbulence can be well described by a 
power-law spectrum in many cases. Focusing on this re- 
gion we model the EGMFs as Kolmogorov turbulence. 
However, these models do not consider magnetic fields 
at spatial scale larger than that where the power-law 
spectrum is valid. The ignorance of the larger scale mag- 
netic fields let us underestimate the deflection angles of 
UHECRs. Since the spectral shape of the magnetic fields 
depends on their generation mechanism, scale and geom- 
etry, we avoid the uncertainty and simply employ Kol- 
mogorov turbulence with the maximum scale so larger 
than the correlation lengths A c and Af, as to reproduce 
the analytical estimation of UHECR deflection angles in 
the last subsection on average. We set A Cjmax = 400 kpc 
and Af jmax = 400 kpc for the maximum scale. 

In any case, once a UHECR is emitted from a source, 
the propagation of the UHECR is calculated by solving 
the equation of motion in each magnetic structure tak- 
ing energy-loss processes by interactions with ambient 
photon fields, i.e., photomeson production and Bethe- 
Heitler pair creation, into account. An event generator 
SOPHIA is adopted to simulate the photomeson produc- 
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Figure 4. Normalized time-delay distributions of protons with 
energies of 10 20 eV produced during propagation in the GMF (ma- 
genta) and extragalactic magnetic structures around a source, i.e., 
a cluster of galaxies when the source is located at the center (red), 
a cluster of galaxies if there is the source at 1 Mpc away from the 
center (green) and a filamentary structure (blue). Note that the 
10 20 eV is the energy of protons at the local frame of the magnetic 
structures. It is the energy of protons at UHECR generation in the 
former three cases, while it is the energy when UHECRs penetrate 
into the Galactic space, i.e., observed energy, in the latter case. 



tion (jMiicke et alJ l200Q) , and the energy- loss rate of pro- 
tons by the Bethe-Heitler pai r creation is estimated by a n 
analytical fitting formula in iCho dorows kTet al.l fl992). 
When the UHECR reaches the boundary of the struc- 
tures, we record the time-delay and deflection angle of 
the UHECR, and obtain these distributions from the 
propagation results of many UHECRs. 

The time-delay and deflection angles of UHECRs 
in the GMF is estimated by usin g a backtrack- 
ing m ethod (e.g.. ITakami fc Sato! [20081: ITakami fe Satol 
20101). For the GMF model, the BS-S model used in 
Takami fc Satol (I20T0I) whi ch was originally proposed by 
Alvarez-Muniz et all (|2002t) . is adopted. Since the GMF 
is distributed similarly to the Galactic arm, the time- 
delay of UHECRs depends on the arrival directions of 
UHECRs. So, we adopt a value averaged over the whole 
sky as a typical value of time spread. All the energy-loss 
processes can be neglected because the propagation path 
length of protons is much smaller than their attenuation 
length at energies we are interested in. 

3.2. Time- delay and deflections 

Figure @] shows the distribution of the time-delay of 
protons with the energy of 10 20 eV in the different 
magnetic structures. o~{E) ~ id{E) is confirmed in 
all the cases. The deviation of the time profile of a 
UHECR burst produced by these magnetic structures 
are a G (E) ~ 10 1 ' 5 yr, a m (E) ~ 10 25 yr, a cl (E) ~ 10 5 
yr and a cc (E) ~ 10 5 8 yr at E = 10 20 eV, respectively. 
Here <tg(E) and <Jfu(E) are the time spread produced 
by the GMF and the EGMF in a filamentary structure, 
respectively. Both a cc (E) and o~ c \{E) are time spreads 
produced by magnetic field in a cluster of galaxies, but 
the former is the case when a UHECR source is located 
at the center of the cluster, while the latter is the case 
when a UHECR source is located at 1 Mpc away from 
the center of the cluster. os\(E), cr c i{E) and a cc (E) are 
proportional to E~ 2 because the turbulent fields are as- 
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Figure 5. Normalized deflection-angle distributions of protons 
with energies of 10 20 eV produced during propagation in the differ- 
ent magnetized regions corresponding to the structures considered 
in Figure [4] The definitions of the labels and the energy of f 2<) 
eV are the same as in Figure [4] 

sumed, which were able to be confirmed by the numerical 
simulations. 

Here, let us check the validity of ip ~ 5°. Figured] 
shows the distribution of the deflection angles of pro- 
tons with the energy of 10 20 eV in the different mag- 
netic structures. The deflection of the protons by the 
GMF and magnetic field in a filamentary structure sat- 
isfy V , ^5°. On the other hand, a cluster of galaxies 
produces the deflection of protons larger than 5°. How- 
ever, in this case, the sky region where UHECRs from the 
source are occupied is determined not by the deflection 
angles but by a viewing angle of the magnetic structure; 
i.e., tp < 5° is satisfied if D > 30 Mpc. Thus, only the 
Virgo cluster (~ 16 Mpc) does not satisfy this require- 
ment. If the Virgo cluster does not have UHECR sources, 
tp ~ 5° is still justified. Otherwise, it is sufficient to con- 
tinue discussions that a larger tp < 10° is adopted. In the 
latter case, the validity of equation (J7J is a bit changed 
as shown in equation ©. We use tp = 5° to simply com- 
pare the case of clusters of galaxies with the other cases 
in this paper. 

4. IMPLICATIONS FOR PROPERTIES OF UHECR 
BURSTS 

For UHECRs observed at energies E, while the time 
spread by the GMF should be estimated at the observed 
energy E, that by structured EGMFs embedding their 
sources should be considered by using E g (E,D), which 
is the energy of UHECRs at generation, because UHE- 
CRs lose their energies during propagation in intcrgalac- 
tic space. This is specially important when E is higher 
than the threshold of the photomeson production. As- 
suming that UHECR sources are uniformly distributed 
in local Universe, we introduce characteristic time spread 
in the structured EGMFs as 



A, 



dD AnD 2 (j yi {Eg{E, D),D), 



47TA„ax (E) 

(15) 

where x = fil, cc, cl. Here E g (E,D) is calculated 
by the backtracking method used for the calculations 
of D maK (E). We adopt the characteristic time profile 



Takami & Murase 



spread t x (E) to constrain p s below. Since E is lower 
than E g (E,D) for any E and D because of energy- loss, 
Cx(-Emax) is generally smaller than t x (E). Therefore, 
Cx(-Emax) instead of t x (E) gives a more conservative con- 
straint. 

The minimum of the characteristic time spread T m i n (E) 
is estimated for the three cases of source locations. If 
the sources are located in filamentary structures, the 
characteristic time spread is estimated as r min (E) « 
max(crG(£ , ),T m (£')) w T m (E), where r m (E) ~ 57 yr at 
E = 10 20 eV. For the remaining two cases where UHECR 
sources are embedded in clusters of galaxies, the time 
spread produced by the clusters is much longer than that 
produced by the GMF, i.e., r min (E) w t cc (E) - 9 x 10 4 
yr and t c1 (E) ~ 2 x 10 4 yr at E = 10 20 eV, respec- 
tively. These values constrain p s from an upper side in 
the bursting case. 

As mentioned before, the EGMF in voids may mainly 
contribute to the total time spread of UHECR bursts. 
Following equations ([3]) and and ip = 5°, the charac- 
teristic time spread by the void EGMF is limited as 



T v {E)^-a v (E,D max {E)) 

< 1 x io 6 e 20 - 2 ( £m^l 

20 \ 75 Mpc 



<im) 



Here, t v (E) dominantly determines T max (_E) in all the 
cases, when t v max{E) is larger than Tfu(E), t c \{E), and 

Tcc(E). 

Then, we can constrain the rate of UHECR bursts fol- 
lowing equation (J7J) for the bursting case, and otherwise 
p s is constrained by equation (fTU)) . into which r max (i?) 
is substituted instead of t(E). Emissivity of UHECRs 
at 10 19 eV required to reproduce the observed UHECR 
flux i s Com ~ 10 44 erg Mpc~ 3 yr" 1 dWaxman fc Bahcalll 
[19991 IBerezinskv eTaLl [2001 IMurase fc Takamil I2009D 
This emissivity can be converted into the requirement 
of cosmic ray energy input per burst by using p a , 
= E 2 (dN^/dE) = Ccr/Ps- Figure 1 summa- 
rizes constraints on p s and of transient UHECR 
sources as functions of n s (E). The black solid line 
represents a constraint of p s derived from T v max (£'). 
The three inclined lines show constraints on p s or 
<?CR- The vertical black dashed line is n s (E) = 3 x 
10- 4 V> 5 ~ 2 (Anax(£)/75Mpc)- 3 Mpc" 3 for E = 10 20 eV 
(see equation (jSjj), which divide the bursting case from 
the other. 

5. DISCUSSION 

Although we focused on transient UHECR sources in 
this work, one should keep in mind that the information 
on n s {E) itself is relevant, whether sources are steady 
or transient. Since actual n s (E) must be smaller than 
nu at any energy E (unless multiple images caused by 
sufficiently strong EGMFs lead to difficulty in determi- 
nation of plausible n s (E)), comparing n s (E) with the 
number density of know astrophysical sources such as FR 
I/II galaxies should be useful, as have been discusse d for 
steady UHECR sources (e.g.. iTakami fc Satdl2009l ). A 
recent constraint on the apparent source number densit 
is - 10" 4 Mpc" 3 above ~ 6 x 10 19 eV (|Takami fc Sat. 
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Figure 6. Possible constraints on the (differential) cosmic-ray 
energy input at 10 19 eV per burst E 2 (dN^/dE) and the rate 
of bursts or flares p s as functions of the apparent source number 
density n s (E). The vertical dashed line distinguishes between the 
bursting case and the other. In the bursting case (the left side of the 
line) these quantities are limited from both sides by equation |[7jl- 
The black solid line shows the left inequality of equation 10, in 
which r max (i?) is dominated by the effective time spread of UHE- 
CRs by the void EGMF t v (E). The three color lines correspond to 
the right inequality of equation (0 for each magnetic environment 
around UHECR sources, i.e., filamentary structures (red), the case 
where sources exist at the center of clusters of galaxies (blue), and 
the case where sources are located at 1 Mpc away from the center 
(magenta). On the other hand, at the right side of the dashed line , 
these quantities are limited only from one side by equation 1101 . 
Note that E = 10 20 eV, B max = 10 21 eV and if> = 5° are assumed. 
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2009; Cuoco ct al. 2009). This value can be regarded as 
n s (E) at E ~ 6 x 10 19 eV because of the steep spec- 
trum observed a t around this energy w ith the spectral 
index of ~ 4.3 (jAbraham et alj I2010bh . Thus, for in- 
stance, FR II galaxies, whic h has the local numbe r den- 
sity of - 3 x 10~ 8 Mpc~ 3 (jBlandford et al.lH990l ). and 
flat-spectrum radio quasars which correspond to on-axis 
FR II galaxies in the framework of the unific ation sce- 
nario of radio galaxies ([Urrv fc Pa dovani 1995), seem too 
rare as UHECR proton sources. 

In this study, we adopt simple three zone models of 
magnetic fields (the GMF, a structured EGMF around 
sources, and the void EGMF), but modeling of the 
EGMFs have large uncertainty at present. Cosmologi- 
cal structure formation simulations have indicated that 
cosmic magnetic fields have complex web structures fol- 
lowin g matter di stribution (jSigl et al.ll2003t IDolag et al.1 
I2005t iDas et all 12008ft . At present, the volume filling 
fraction of relatively strong (> 10 nG) magnetic fields, 
i.e., filaments and cluster s, highly depends on mo dels and 
methods (see Figure 9 of IKotera fc Olintol ([201 If )). Also, 
the spectrum of the magnetic fields is crucial for the prop- 
agation of UHECRs, though it is not certain. An obser- 
vational constraint of the magnetic stren gth of EGMFs in 
filamentary structures is only < 0.1/nG ([Xu et al.|[2006ft 
and there is no direct observational implication on the di- 
rection and coherent length of the magnetic fields. While 
a turbulent magnetic field with the Kolmogorov spec- 
trum with Bf = 10 nG is simply adopted in this work, 
several numerical simulations have indicated the exis- 
tence of large-scale fi elds in filamentary structures (e.g., 
Briiggen ct al. 2005). Such a coherent component de- 
flects the trajectories of UHECRs more efficiently and 
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Figure 7. Same as Figure[6] but the possible case that a magnetic 
field in a local filament contributes to the propagation of UHE 
protons if UHECR sources are embedded in filamentary structures 
(blue). For reference, a limit when only the GMF contributes to the 
propagation, i.e., the effect of a magnetic structure around sources 
is negligible, is s hown (red), which is reg arded as a conservative 
case discussed in Murase & Takami (2009). The black dashed line 
and solid line are the same as in Figure[6] Note that, for the cases 
where sources are embedded in clusters of galaxies, constraints on 
p B and £q)^ given in Figure [6] are more relevant. 



makes longer time-delay compared to a turbulent com- 
ponent. In this case the time spread of a UHECR burst 
is smaller than the time-delay because particles similarly 
propagate due to the large-scale component. 

The complex web of EGMFs implies that structured 
EGMFs on the way from sources to the Milky Way out- 
side EGMFs embedding the sources could also play a 
significant role in the total deflection and time sperad 
of UHECRs, but these do not affect our constraints de- 
rived in the previous section in the scope of our EGMF 
models. Possible effects of such EGMFs are making ad- 
ditional deflection and time spread of UHECRs. As long 
as the additional deflection is small enough, the time 
spread by the EGMFs embedding the sources provides 
conservative constraints. In order to see the effects, we 
begin from evaluating the probability that UHECRs en- 
counter a magnetic structure on the way to the Earth 
during propagation, which can be calculated from the 
number density, n x , and cross-sectional area, A x , of the 
structure x. The number density, n x ss /x/Kc, where V x 
is the volume of the structure x, is rif « 3 x 1CP 5 Mpc -3 
and n c « 10 -6 Mpc -3 , respectively. The cross-sectional 
areas are A f « 2 2 x 25 = 100 Mpc 2 and A c « 7r x 3 2 = 9n 
Mpc 2 , so the encounter probability is 0.24 and 2 x I0~ 3 
for 75 Mpc propagation in the cases of filaments and 
clusters, respectively. Thus, whereas it is unlikely for 
UHECRs to encounter clusters during their propagation, 
some of UHECRs penetrate into a filamentary structure 
once. But their typical deviation angle of protons with 
energies above 10 20 eV is still smaller than the considered 
value of ip = 5° (see Figure O, so that such intervening 
EGMFs do not affect our conservative results. Stronger 
constraints on p s and Sqr from additional time spread 
can be expected especially if UHECRs must pass the lo- 
cal magnetized structure, as discussed afterwards. 

One might expect the situation where effects of struc- 
tured EGMFs are even more prominent. We could mis- 



perceive the positions of UHECR sources from their ar- 
rival directions if UHECRs propagate selectively along 
the magnetic web structure and/or are strongly scat- 
tered off by the struct ured EGMFs (jKotera fc Lemoinel 
120081 IRvu et aUkOlOf ). When such effects are large in 
the Universe, taking larger values of if) will be more 
realistic. Also, as recently suggested, the structured 
EGMFs may increase the probability that we observe 
UHECRs from transient sources thanks to significant 
time spread longer than the intrinsic burst duration of 
the sources (|Kalli et al.ll201ll ). 

The Milky Way is thought to be located in a dense 
region in local Universe. Although magnetic fields in the 
local structure is highly uncertain, it could be another 
unavoidable EGMF for UHECRs arriving at the Earth. 
A constrained simulation has shown that the Milky Way 
belongs to a filamentary structu re adjacent to the local 
Supercluster (jKlypin et al.l 120031 ). Fl Assuming our fila- 
ment model for the local magnetic structure, this EGMF 
is dominant in T m i n (E) in the case where sources are lo- 
cated in filamentary structures because the time profile 
spread should be estimated as <jf(E) instead of Tf(E). 
Figure [7] shows possible constraints on p s and £gg in 
the case where the local filament is included. Conser- 
vative constraints by th e GMF, which is discussed in 
iMurase fc Takanil (|2009j) . is also shown. For a filamen- 
tary structure between a source and the local environ- 
ment discussed above, the time profile spread can also 
be conservatively estimated by af(E), and therefore the 
same constraints as the case of the local filament are 
applied even if UHECRs pass several filamentary struc- 
tures. In the cases that sources are embedded in clusters 
of galaxies, constraints are unchanged because the time 
spread of UHECRs by EGMFs in the clusters of galaxies 
is much larger. The discussions in Section 01 where only 
EGMFs surrounding the sources are considered, provide 
conservative constraints. If the volume filling fraction of 
structured EGMFs and the nature of local magnetic en- 
vironments are understood, we would be able to obtain 
better constraints. 

As seen above, it is obviously crucial to reduce the 
uncertainty of the EGMFs both theoretically and ob- 
servationally in future to understand the properties of 
transient UHECR sources. In order to get better knowl- 
edge on the EGMFs, Faraday rotation surveys by fu- 
ture experiments such as Square Kilometer Arrajlj are 
useful. Detection of TeV synchrotron pair echo/halo 
emission produced by UHE gamma-ray bursts/flares 
also may allow us to probe structured EGMFs with 
> 10 nG 4M urasd l20il[) . Future numerical simulations 
on structured EGMFs can also give us more insight into 
consequences to transient UHECR source population. 

Throughout this paper we have focused on protons. 
Even though heavy nuclei are considered, the discussions 
in this paper are in principle applicable. A basic differ- 
ence between protons and heavy nuclei is electric charge, 
i.e., nuclei suffer from deflections Z times as large as 

4 A larger local magnetic field with ~ 0.1 p,G may be possi- 
ble if the Milky W ay is located in the edge of the Virgo cluster 
(Blasi ct al. 1999). Such a strong local magnetic field may be rel- 
evant as well as the void EGMF, especially when the strong effec- 
tive m agnetic field is required as in the scenario of high-luminosity 
GRBs (Murase fc Takami|[2009l) . 

5 http://www.skatelescope.org/ 
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protons and the time spread roughly Z 2 times as large 
as protons. If structured EGMFs only in the vicinity of 
UHECR sources affect the propagation of nuclei, the dis- 
cussions are valid because the deviation scale, ip, reflects 
not the deflection angles of nuclei in the magnetic struc- 
tures but a viewing angle of the structures. However, 
the GMF and a magnetic field in the Local Group would 
become seriously important. For the GMF, since the de- 
flections of nuclei are too large at energies less than 10 20 
eV, we should focus on nuclei with energies much higher 
than 10 20 eV. The deflection angles of nuclei could be 
the order of tens deg ree above a few times 10 20 eV (e.g., 
iGiacinti et "ail (|201Cfl ). Thus, even considering the GMF, 
the discussions can in principle possible with larger val- 
ues of %/} (e.g., ip ~ 30°). It is very uncertain that such 
discussions are possible even when an EGMF in local 
Group is taken into account because of more uncertainty. 
Future studies on effects of the EGMF in local Group will 
be necessary. 

Although we have focused on an extragalactic 
component of U HECRs, possible Galactic sources, 
includin g GRBs (iLevinson fc Eichlerl 119931 ). hyper- 
novae (IBudnik et al I 120071 ). fast-rotating neutron 
stars (Bl asi et al. 2000) and magnetars, also could con- 
tribute to observed UHECRs. In general, Galactic source 
scenarios have difficulty to reproduce a nisotropy in the 
arrival distribution of UHECRs (e.g ., iPohl fc Eichlerl 
120111 ). Recently iCalvez et all (|2010| ) suggested that 
Galactic sources would cause the gradual change of 
UHECR composition reported by the PAO and therefore 
also contribute to the UHECR flux around the ankle 
without inconsistency with observed anisotropy. How- 
ever, this scenario does not affect our discussions, since 
it still requires an extragalactic component of UHECRs 
above 3 x 10 19 eV, which we are interested in. 

Although we have considered how we can reveal 
transient sources with observations of charged cos- 
mic rays, for transient sources, the multi-messenger 
approach is definitely rele vant, i.e., high-energ y neu- 
trino observations (e.g., iWaxman fc Bahcall I 119971 : 
Rachen fc Mcszaros 1998), Ge V-TeV gamma-r a ys (e.g., 
Atovan fc Dermer I 1200% Dermer et all l2009t iMurasd 
2011J), and UHE photons (|Murasd 120011 are helpful for 
identifying the sources. 

6. SUMMARY 

We studied the propagation of UHECRs with analyti- 
cal models of structured EGMFs, and demonstrated the 
effects of the structured fields on revealing the properties 
of transient UHECR sources by UHECR experiments. 
While the void EGMF seems difficult to probe, struc- 
tured EGMFs may be measured by future radio obser- 
vations or revealed by dedicated numerical simulations. 
Then, UHECR experiments with large exposures may 
allow us to constrain transient UHECR source popula- 
tion by comparing the derived properties, i.e., p s and 
^CRi with those of known astrophysical transients such 
as GRB s, AGN flares and magne t ar generation (see ta- 
ble 4 of iMurase fc Takamil (|2009t ) for generation rates). 
We demonstrated that the energy- dependence of the ap- 
parent source number density n s (E) is a crucial hint 
of transient sources, where, as we argued, observations 
above 10 20 eV are desirable due to the possible multiple- 



burst contamination at lower energies. Once this tran- 
sient feature is identified, the rate of UHECR bursts p 8 
and energy input per burst £q^ can be limited, given 
good knowledge on the structured EGMFs. 

Finally, based on the discussions in this paper, we sug- 
gest a strategy to identify transient UHECR source pop- 
ulation by future UHECR experiments with large expo- 
sures: 



1. to estimate n s {E) from high event statistics and to 
see whether UHECR sources are transient from the 
energy dependence of n s (E) , 

2. to calculate the time profile spread of UHECRs 
from bursting or flaring sources, given good knowl- 
edge on the structured EGMFs from observations 
and/or simulations 



3. to constrain p s and £(?r from the above pieces of 
information 



4. to discuss implications for known astrophysical 
candidates of UHECR sources via comparison with 
the properties of these transients. 



We thank to P. L. Biermann, C. D. Dermer, and S. In- 
oue for useful comments and discussions. We also grate- 
ful to the anonymous referee. The work of K.M. is sup- 
ported by a Grant-in- Aid from Japan Society for the Pro- 
motion of Science (JSPS) and the Center of Cosmology 
and AstroParticle Physics (CCAPP). 
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